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HIGHLIGHTS 


•  A  model  for  a  water-cooled  PEM  fuel  cell  stack  was  developed. 

•  A  time-varying  PID  controller  was  implemented  on  this  model. 

•  The  controller  was  tested  on  an  experimental  power  profile  from  a  3-bedroom  house. 
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In  this  paper,  a  time-varying  proportional-integral  (PI)  controller  is  designed  for  controlling  the  tem¬ 
perature  of  a  water-cooled  5  kW  hydrogen  fuel  cell  stack.  This  controller  is  designed  using  a  mathe¬ 
matical  model  for  the  stack  temperature,  which  is  derived  using  basic  chemical  engineering  material  and 
energy  balances.  The  controller  affects  the  stack  temperature  by  changing  the  flow  rate  of  cooling  water 
that  passes  across  the  stack.  The  model  is  then  analyzed  using  a  number  of  power  demand  profiles  to 
determine  the  effectiveness  of  the  controller.  The  results  show  that  a  time-varying  PI  controller  is 
adequate  for  maintaining  the  stack  temperature  within  5  K  of  the  target  point. 

©  2014  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  the  search  for  alternative  fuel  sources,  hydrogen  has  become  a 
popular  option  primarily  due  to  the  advantages  offered  by 
hydrogen  fuel  cells  over  the  modern  internal  combustion  engine 
[1].  Fuel  cells  have  multiple  advantages  over  conventional  internal 
combustion  engines  such  as  improved  efficiency,  lack  of  moving 
parts,  reduced  noise  during  operation  and  the  ability  to  store  en¬ 
ergy  directly  in  batteries  without  any  additional  steps.  Although 
fuel  cells  are  already  readily  available  for  use  as  backup  generators 
and  other  stationary  power  sources,  the  idea  of  using  a  fuel  cell  for  a 
dynamic  power  demand  is  relatively  new  and  the  development  of 
such  a  system  is  still  being  researched.  One  of  the  biggest  concerns 
in  this  situation  is  the  ability  to  keep  the  fuel  cell  stack  at  optimal 
conditions  to  maintain  the  highest  efficiency  in  the  face  of  fluctu¬ 
ating  power  demand.  In  the  case  of  a  constant  power  demand,  it  is 
much  easier  to  maintain  the  stack  at  optimal  operating  conditions, 
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but  a  dynamic  power  demand  requires  some  reaction  to  the 
changes  so  that  the  efficiency  of  the  stack  can  be  maintained.  In  this 
paper,  we  focus  on  the  dynamic  analysis  and  control  of  a  water 
cooled  fuel  cell  stack,  which  is  being  used  as  a  backup  generator  in  a 
stationary  environment  (e.g.  a  residential  house). 

While  there  have  been  a  large  number  of  papers  in  the  area  of 
modeling  fuel  cells  for  different  applications,  very  few  of  these 
papers  have  focused  on  the  issue  of  temperature  dynamics  of  the 
fuel  cell  and  its  effect  on  performance.  Some  of  the  earliest  research 
into  the  thermal  characteristics  of  a  PEM  fuel  cell  was  actually 
experimental.  Amphlett  et.  al.  [2]  conducted  numerous  physical 
experiments  on  a  working  fuel  cell  stack  and  determined  properties 
such  as  the  effective  mass  and  surface  area  of  the  fuel  cell  stack. 
These  values  compared  well  with  estimates  based  on  bulk  stack 
measurements.  The  fuel  cell  stack  was  also  run  at  various  experi¬ 
mental  conditions  to  determine  the  stack  response.  These  results 
were  then  used  to  experimentally  determine  further  characteristics 
of  the  stack  including  a  prediction  for  the  cell  voltage,  the  heat 
transfer  coefficient,  and  the  thermal  capacity  of  the  stack.  A  tran¬ 
sient  model  was  also  constructed  to  determine  the  stack  temper¬ 
ature  with  time,  but  relies  on  a  large  number  of  inputs  that  are,  at 
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this  point,  experimentally  determined.  Grasser  and  Rufer  [6] 
developed  a  thermal  model  for  liquid  cooling  of  a  PEM  fuel  cell  in 
their  dynamic  model.  With  their  experimental  results,  they  were 
able  to  design  the  thermal  model  using  simple  energy  balances  to 
express  the  stack  temperature  in  a  first-order  differential  equation, 
similar  to  how  Amphlett  et  al.  performed  their  research. 

A  variety  of  theoretical  models  have  been  developed  to  study 
and  model  PEM  fuel  cells.  Pukrushpan  et.  al.  [12]  developed  a 
nonlinear  dynamic  model  that  is  suitable  for  control  study.  Xue  et. 
al.  14]  developed  a  dynamic  lumped-parameter  model  of  a  PEM 
fuel  cell  and  discovered  that  hydrogen  diffusion  through  the  elec¬ 
trode  affects  the  stack  temperature.  Pathapati  et.  al.  [11]  further 
developed  the  model  set  forth  by  Xue  et.  al.  [14]  specifically 
focusing  on  the  dynamics  and  mechanisms  that  are  of  particular 
importance  for  automotive  purposes.  Chu  et.  al.  [3]  also  developed 
a  lumped  parameter  model,  and  investigated  the  changes  in  stack 
temperature  for  various  inputs,  but  control  of  the  stack  tempera¬ 
ture  was  not  in  the  scope  of  their  research.  Meyer  and  Yao  [9] 
presented  a  Multiple  Input  Single  Output  (MISO)  cooling  system 
model  using  linearization  and  ^-synthesis,  but  this  model  is  again 
rather  complicated  and  integrated  into  an  overall  dynamic  model 
for  the  stack.  Other  research  has  been  performed  to  analyze  how 
the  fuel  cell  stack  with  perform  under  adverse  conditions,  such  as 
starting  the  stack  in  freezing  temperatures  [13]. 

While  various  theoretical  models  have  been  developed  to  study 
the  PEM  fuel  cell  stack,  they  tend  to  have  thermal  models  devel¬ 
oped  for  a  fuel  cell  stack  open  to  the  atmosphere.  This  air-cooled 
stack  is  of  interest  in  the  laboratory,  but  it  does  not  realistically 
model  what  the  temperature  profile  of  the  stack  would  begin  to 
look  like  when  used  in  an  automobile.  When  in  operation,  the  at¬ 
mosphere  under  the  hood  of  a  vehicle  can  get  much  hotter  than  the 
surrounding  environment  outside  of  the  vehicle,  making  the  air  in 
the  engine  compartment  unsuitable  for  adequate  cooling  of  the  fuel 
cell  stack.  If  an  air-cooled  system  is  not  possible,  the  only  other 
ways  to  reliably  affect  the  stack  temperature  are  to  modify  the  inlet 
conditions  or  to  use  a  liquid  cooling  system  around  the  outside  of 
the  stack.  Since  modifications  of  the  inlet  conditions  will  obviously 
have  an  adverse  effect  on  the  desired  power  output,  a  liquid  cooling 
system  is  the  only  other  feasible  option.  Significantly  less  theoret¬ 
ical  research  has  been  done  on  liquid  cooling  of  a  fuel  cell  stack. 
Recently,  Fang  et.  al.  [5]  developed  a  dynamic  model  using  a  liquid 
cooling  thermal  model,  but  control  of  the  stack  temperature  was 
once  again  not  in  the  scope  of  their  research. 

In  this  research,  a  theoretical  model  is  developed  to  model  the 
thermal  characteristics  of  a  hydrogen  fuel  cell  stack,  as  well  as  its 
response  to  changes  in  parameters  that  can  be  used  in  temperature 
control.  A  schematic  of  the  system  under  consideration  is  shown  in 
Fig.  1. 

The  specific  objectives  are  (1)  to  develop  a  time-varying  pro¬ 
portional-integral  (PI)  controller  for  cooling  water  flow  rate  across  a 
jacketed  PEM  fuel  cell  stack  with  the  intent  of  controlling  the  stack 
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temperature  and  (2)  test  the  controller  using  a  variety  of  power 
profiles. 

Fuel  cell  stacks  currently  offer  one  of  the  best,  cleanest  alterna¬ 
tives  to  the  fossil  fuels  used  in  modern  engines,  thereby  reducing  the 
amount  of  harmful  greenhouse  gases  released  into  the  atmosphere. 
Improvement  in  the  design  and  efficiency  of  fuel  cell  stacks  will  also 
make  them  more  affordable  for  stationary  power  applications. 

2.  Development  of  thermal  model 

A  fundamental  premise  of  the  approach  taken  here  is  that  the 
electrical  characteristics  of  a  fuel  cell  change  much  more  rapidly 
than  the  temperature.  A  time-scale  analysis  of  this  system  indicated 
that  the  electrical  time-scale  was  of  the  order  of  milliseconds,  the 
mechanical  flow  process  time-scale  was  of  the  order  of  seconds  and 
the  thermal  time-scale  was  of  the  order  of  minutes  [8  ].  Thus  we  can 
examine  the  dynamics  of  the  stack  temperature  assuming  that  the 
voltage  adjusts  nearly  instantaneously  when  current  draw  changes. 
This  does  not  mean  that  the  voltage-current  relationship  does  not 
depend  on  temperature,  just  that  the  changes  occur  so  rapidly  that 
this  part  of  the  overall  model  can  be  represented  by  an  algebraic 
model. 

If  the  fuel  cell  stack  and  its  contents  are  chosen  as  the  system, 
the  energy  balance  can  be  written  as: 

=  E  ^  inW/  in  -  E  ^outHi.out  ~Q.~W  ( 1 ) 

i  i 

The  material  balance  for  any  species  i  is  the  following: 


=  Ain out +  «'!?  (2) 

Even  though  there  are  separate  reactions  occurring  at  the  anode 
and  cathode,  we  combine  them  into  an  overall  water  formation 
reaction.  The  accumulation  of  the  species  i  is  very  small  compared 
to  the  mass  of  the  stack.  It  is  therefore  safe  to  assume  steady-state 
for  all  gaseous  and  liquid  species.  This  assumption  allows  the 
accumulation  term  in  the  material  balance  to  be  set  equal  to  zero, 
resulting  in: 

M/,out  <m  t^i, in  +  vi%  (3) 

The  material  and  energy  balances  may  now  be  combined  to 
give: 

^  =  E  N-  in  (Hi, in  -  Hi, out)  -  AHrxn(r0Ut)|  -Q  -W  (4) 


where  the  heat  of  reaction  for  the  overall  reaction,  the  formation  of 
water,  will  be  evaluated  at  the  actual  cell  outlet  temperature.  Note 
that  since  the  NI  0Ut  term  can  be  written  in  terms  of  the  Ni  in  term 
from  Eq.  (3),  the  combined  equation  presented  in  Eq.  (4)  is  a 
function  of  Ni  in  only. 

If  the  hydrogen  for  the  fuel  cell  is  generated  by  the  reforming  of  a 
hydrocarbon  then  the  anode  inlet  flow  has  hydrogen  (H2),  carbon 
dioxide  (C02),  carbon  monoxide  (CO)  and  water  (H20).  The  cathode 
inlet  flow  is  made  up  of  oxygen  (02),  nitrogen  (N2),  and  water  (H20). 
All  required  species  flows  can  be  related  to  the  inlet  hydrogen  flow. 
The  mole  fraction  of  hydrogen  in  the  inlet  is  given  by: 


y  h2 


N 


H2,in 


^H2,in  +  Mother, in 


(5) 


Then  all  of  the  other  gas  flows  (carbon  dioxide,  carbon  mon¬ 
oxide,  and  water)  can  then  be  expressed  by 


Fig.  1.  Basic  schematic  of  fuel  cell  system. 
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Nother.in  =  (^^)NH2,in  (6) 

The  cathode  inlets  can  similarly  be  expressed  as  a  function  of  the 
hydrogen  inlet  flow. 


r02,in  =  (1  +Xs)NH2,in 

(7) 

rN2,in  -  §(1  +Xs)JVH2,in 

(8) 

W  =/s(l+Xs)NH2,in 

0) 

Voell  =/(/,Tst)  (15) 

where /is  a  functional  form  representing  the  polarization  curve.  In 
a  fuel  cell  stack,  the  polarization  curve  captures  the  reduction  in 
voltage  with  increasing  current  density  due  to  ohmic  and  transport 
resistances  and  is  a  function  of  temperature,  pressure  and  stack 
material  properties.  For  maximum  efficiency,  it  is  desirable  to 
operate  at  a  point  in  the  polarization  curve  where  the  power  pro¬ 
duced  is  maximum.  Thus,  a  change  in  stack  temperature  changes 
the  operating  point  on  the  polarization  curve.  Combining  all  of 
these  equations,  the  overall  energy  and  material  balance  now  be¬ 
comes  the  following. 


where  fs  represents  the  humidity  ratio  of  the  entering  air.  The 
previous  equations  also  take  into  account  any  excess  of  oxygen  in 
the  system  and  a  cathode  stream  that  is  saturated  with  water.  The 
saturation  calculation  may  be  solved  or  estimated  as  necessary.  This 
information  can  now  be  used  to  simplify  the  inlet  flow  terms  of  the 
combined  material  and  energy  balance. 

E  E"  (Hi, in  -  Hi, out)  =  |  (T-jn  -  Tstack)Cp,m  ( 1 0) 

i 

where 


-p,m 


Cp,H2  +  (  ytfH2)Cp,CQ,c°2,H20  +  ^ 

+  nj(  1  "T  Xs)Cp;N2  +/S(1  +  ^s)Cp,H20 


(11) 


and  X  represents  the  fractional  conversion  of  hydrogen  in  the  cell 
and  Cp  m  is  an  effective  heat  capacity  for  the  entering  gases.  The 
inlet  flow  term  can  further  be  simplified  by  setting  a  constant 
a  =  Cp,m/X.  The  inlet  flow  rate  terms  can  now  be  expressed  as: 

E^4n(HUn  -  tffout)  =  afOin  -  Tst)  (12) 

i 

The  change  in  extensive  energy  term  in  the  combined  energy 
and  material  balance  can  also  be  simplified.  It  is  known  that  the 
extensive  energy  of  the  system  is  equal  to  the  sum  of  the  extensive 
energy  of  the  stack,  and  the  extensive  energy  of  the  contents  of  the 
stack.  It  can  again  be  assumed  that  the  mass  of  the  stack  makes  any 
change  in  the  contents  relatively  insignificant  when  it  comes  to  the 
stack  temperature.  The  differential  term  can  therefore  be  approxi¬ 
mated  as  the  change  in  the  extensive  energy  of  the  stack  [2]. 


d U 
dt 


ttlstCp^st 


dTst 

dt 


(13) 


where  ms t  and  CPiSt  represent  the  mass  and  specific  heat  of  the  fuel 
cell  stack.  The  work  term  of  the  combine  energy  and  material 
balance  must  be  defined.  Since  the  work  done  by  the  stack  is 
electrical  power,  it  is  only  logical  to  express  it  as  such.  The  power 
output  for  each  fuel  cell  is  the  product  of  current  draw  and  voltage. 
The  voltage  is  a  function  of  current  draw,  current  density,  and  stack 
temperature.  The  total  power  output  for  the  system  is  therefore  the 
power  output  from  each  cell  multiplied  by  the  number  of  fuel  cells 
in  the  stack. 


dT  £ 

mstCp,sf-a^=J(Tin-rst)Cp,m-AHrXn(Tout)?-Q-nceiis/Vcell(/,Tst) 

(16) 

The  heat  of  reaction  term  shown  in  Eq.  (16)  must  be  calculated 
at  the  actual  operating  conditions  for  the  fuel  cell  stack.  This  means 
that  it  must  be  evaluated  at  the  stack  operating  temperature  and  for 
the  appropriate  state  (liquid  or  vapor)  of  the  reaction  products. 
Because  we  can  be  in  a  situation  where  some  of  the  water  formed 
can  be  in  liquid  state,  the  actual  heat  of  reaction  can  be  calculated 
from 


AHrxn  —  /l^^rxn.l  +  0  — /l)^^rxn,v  —  AHrXn,v  — /AHVap  (17) 


where /i  is  the  fraction  of  the  water  formed  that  exits  as  a  liquid,  the 
subscripts  on  the  heat  of  reaction  terms  refer  to  the  values  calcu¬ 
lated  for  the  water  produced  as  either  a  vapor  or  a  liquid,  and  the 
term  AHvap  is  the  heat  of  vaporization  of  water.  To  justify  neglecting 
the  last  term,  we  should  realize  that  the  ratio  of  the  heat  of 
vaporization  to  the  heat  of  reaction  of  the  vapor  is  given  by: 


AiTvap  44 
AiTrxn,v  241.8 


0.182 


(18) 


This  if  even  as  much  as  10%  of  the  water  exits  as  a  liquid  (this 
would  be  a  very  high  value),  the  possible  error  in  using  only  the 
heat  of  reaction  to  form  water  vapor  (and  neglecting  the  /iAHvap 
term  in  Eq.  (17))  is  less  than  2%.  A  more  detailed  model  could 
incorporate  the  phase  equilibrium  calculation  that  would  allow  for 
the  calculation  of  the  fraction  that  exits  as  a  liquid.  We  hope  to 
address  this  point  in  a  future  publication. 

Eq.  (16)  can  be  further  simplified,  however,  since  the  extent  of 
reaction  can  be  related  to  current  draw. 


I  =  fil  (19) 

where  denotes  the  conversion  from  material  flow  to  current.  This 
can  be  calculated  using  Faradays  constant,  F,  and  the  number  of 
electrons  involved  in  the  reaction,  ne  .  When  substituted  into  the 
material  and  energy  balance  and  combining  all  of  the  constant 
parameters  into  singular  constants,  the  inlet  and  outlet  flow  terms 
can  be  set  as  a  function  of  current  draw. 

mstCp.st^f  =  C,I(Tin  -  Tst)  -  C2I  -  Q  «  nceUsIVceU  (20) 

The  constants  Cj  and  C2  can  be  calculated  as  follows 


w  =  nCei|/Vceii  (14)  d  =  TwiL 

xne-  r 

where  the  cell  voltage  VCeii  is  a  function  of  the  current  density  (z  =  // 

Acen)  and  the  stack  temperature  (Tst)  that  is  represented  by  the  r  _  AHrxn 
polarization  curve.  Thus,  we  can  express  2  ~  ne-F 


(21) 

(22) 
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Just  like  the  fuel  cell  stack  itself,  the  cooling  water  used  to 
regulate  the  stack  temperature  has  its  own  dynamics.  However,  the 
cooling  water  temperature  is  also  dependent  on  position.  If  the 
cooling  water  jacket  is  divided  into  n  zones,  where  the  temperature 
within  each  zone,  i,  is  assumed  to  be  uniform,  the  energy  balance 
for  each  zone  can  be  defined  as: 


dr 

dt 


the  w 

pvT 


(Ti- 1 


where 


(23) 


Qi  =  VMTst  -  Tt)  (24) 

This  is  analogous  to  considering  a  plug-flow  reactor  as  a  number 
of  CSTRs  in  series.  This  simplification  results  in  a  model  for  the 
cooling  water  temperature  dynamics  that  is  a  series  of  ordinary 
differential  equations.  This  simplification  essentially  implies  that 
the  stack  itself  is  at  a  uniform  temperature  and  with  this  assump¬ 
tion,  there  is  no  effect  of  co-current  or  counter-current  flow.  While 
this  “lumped  model”  approach  might  not  be  accurate  for  a  very 
large  stack  (e.g.  one  that  supplies  100  kW  of  power),  for  a  small 
stack  that  produces  up  to  5  kW  of  power,  this  is  a  reasonable 
approach  as  shown  from  the  experimental  data  of  Khan  and  Iqbal 
[7]. 

The  heat  transfer  term  for  the  fuel  cell  stack  can  be  shown  to  be: 


d  =  -  E  =  -  E  UA(T*  -  5)  (25) 

i  =  i  i  =  i 

This  results  in  the  following  stack  temperature  dynamics: 


mstCPst 


dTst 

dt 


Xne-F 


/(Tin  -  Tst)  - 


AHn 


ne-F 


l-UA^Tst-T,) 

1  =  1 


^cell^cell 


(26) 


A  proportional-integral  (PI)  controller  is  modeled  in  state  space 
as: 


rhCw  —  thc  w,ss  +  kc(Tsp  —  Tst)  +—  J  (Tsp  —  Tst)dt  (27) 

o 

where  the  first  term  rhcw, ss  is  the  flow  rate  of  cooling  water  at 
steady-state,  the  second  term  represents  the  proportional  gain 
term  in  the  control  law  and  the  third  term  represents  the  integral 
term  in  the  control  law. 


3.  Simulation  studies 

The  model  developed  earlier  was  simulated  in  MATLAB.  The 
model  parameters  used  are  for  a  Ballard  Mark  V  5  kW  Fuel  Cell 
stack  and  are  given  in  Table  1  [7].  It  is  necessary  to  parameterize  the 
voltage  versus  current  density  relationship  in  16  depending  on 
thermodynamic  state  of  the  system.  This  polarization  curve  is 
parameterized  using  the  approach  proposed  by  Pukrushpan  et  al. 
[12].  A  MATLAB  program  was  written  to  determine  the  voltage  as  a 
function  of  current,  gaseous  partial  pressures  and  the  temperature 
at  each  instance.  More  details  can  be  found  in  Ref.  [10]. 

Before  performing  the  calculations,  it  is  first  necessary  to 
determine  the  number  of  zones  that  the  cooling  water  channel 
should  be  broken  into.  The  number  of  zones  was  determined  by 
performing  the  calculations  using  various  numbers  of  zones  and 
determining  at  what  point  the  temperature  of  the  final  zone  was 


Table  1 

Parameters  from  a  Ballard  Mark  V  5  kW  fuel  cell  stack. 


Cooling  water  specific  heat,  CPcw 

4.1 796  Jg"1  K  1 

Cooling  water  density,  pcw 

1000  kg  irr3 

Inlet  temperature  of  cooling  water,  Tin 

300  K 

Heat  transfer  parameter  for  stack,  UAst 

35.55  +  0.025/  W  K1 

Number  of  cells,  nceii 

35 

Stack  length,  I 

0.38  m 

Stack  width,  W 

0.21  m 

Cooling  water  channel  thickness,  h 

0.02  m 

Product  of  stack  mass  and  specific  heat,  mCst 

35  kj  K-1 

Lumped  parameter  CPgas  /Xne  F 

1.544  WK"1  A-1 

Lumped  parameter  AHrxn/ne  F 

626.4  W  A-1 

Active  fuel  cell  area,  Acen 

232  cm2 

consistent.  There  is  practically  no  change  in  the  temperature  of  the 
final  zone  when  more  than  15  zones  are  considered,  so  15  zones 
were  used  in  this  simulation  study. 

The  first  term  in  the  PI  controller  design  equation  is  the  cooling 
water  flow  rate  needed  to  keep  the  stack  at  the  optimal  tempera¬ 
ture.  The  literature  shows  that  the  target  operating  temperature  for 
a  PEM  fuel  cell  is  about  350  K  [1  ].  The  difficulty  arises,  however,  in 
considering  the  fact  that  the  needed  cooling  water  flow  rate  is 
dependent  on  the  current  draw;  greater  current  draws  will  result  in 
greater  power  output,  but  also  greater  heat  produced.  This  addi¬ 
tional  heat  must  be  removed  from  the  stack  to  maintain  its  tem¬ 
perature  at  the  target  value. 

The  first  task  is  therefore  to  determine  the  relationship  between 
the  current  draw  and  the  cooling  water  flow  rate  to  maintain  the 
stack  temperature  at  350  K.  This  was  done  by  considering  the  same 
model  used  for  calculating  the  stack  dynamics  and  setting  the  stack 
temperature  as  constant.  It  is  then  possible  to  perform  the  calcu¬ 
lations  for  a  series  of  current  draws  and  determine  the  necessary 
cooling  water  flow  rate.  The  needed  cooling  water  flow  rates  for  an 
inlet  cooling  water  temperature  of  300  K  are  plotted  in  Fig.  2. 

A  time-varying  PI  controller  was  designed  with  the  goal  of 
minimizing  the  change  in  stack  temperature  from  the  optimal 
point.  In  quantitative  terms,  the  goal  was  to  develop  a  controller 
that  would  keep  the  stack  temperature  within  5  K  of  the  target 
point.  The  model  was  given  a  number  of  current  draw  profiles  to 
test  the  controllers  ability  to  react. 

The  parameters  used  to  define  the  PI  controller  are  given  in 
Table  2.  These  parameters  were  determined  by  first  considering 
their  proper  order  of  magnitude  such  that  the  first  term  was  the 
major  contributor  to  the  steady-state  value  of  the  cooling  water 


Fig.  2.  Cooling  water  flow  rate  to  achieve  stack  temperature  of  350  K. 
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Table  2 

PI  controller  parameters. 


Set  point,  Tsp 

350  K 

Steady  state  water  flow  rate,  mcw, ss 

2-36  g  s"1 

Controller  gain,  kc 

-0.03  g  s— 1  K1 

/<  25  A 

-0.02  g  s— 1  K1 

/>  25  A 

Integral  time  constant,  ti 

75  s 

/<  25  A 

125  s 

/>  25  A 

flow  rate.  The  constants  were  then  adjusted  using  a  guess-and- 
check  method  until  the  results  of  the  model  had  been  optimized. 
The  values  for  these  parameters  were  then  used  to  attempt  to 
control  the  stack  temperature  for  other  current  draw  profiles. 

It  was  determined  that  a  different  set  of  controller  constants  was 
more  effective  for  low  current  draws  where  the  stack  is  operating  in 
a  non-linear  portion  of  the  polarization  curve.  This  is  due  to  a 
combination  of  the  polarization  curve  being  non-linear  and  the 
small  changes  needed  in  the  cooling  water  flow  rates  compared  to 
the  changes  for  larger  current  draws.  The  needed  cooling  water 
flow  rate  for  current  draws  below  20  A  is  very  small  compared  to 
the  flow  rate  needed  for  larger  current  draws,  meaning  that  when 
the  current  draw  changes,  the  controller  makes  much  smaller 
changes  to  the  flow  rate  than  it  does  when  the  current  draw  is 


greater.  The  controller  constants  are  therefore  adjusted  so  that  the 
second  and  third  terms  are  more  effective  at  adjusting  the  steady- 
state  flow  rate. 

Fig.  3  shows  a  linear  current  draw  (the  simplest  current  draw 
profile  considered)  and  the  effect  of  the  PI  controller  on  the  stack 
temperature  for  this  current  draw  profile  for  a  constant  cooling 
water  flow  rate  of  5  g  s-1,  for  comparison.  Fig.  3  shows  that  the 
controller  is  extremely  effective  at  maintaining  the  stack  temper¬ 
ature  under  such  simple  demands,  with  virtually  no  change  in  the 
stack  temperature  from  the  optimal  value.  The  plots  of  the  corre¬ 
sponding  cooling  water  flow  rate  and  temperature  are  shown  Fig.  4. 

Further  analysis  of  the  results  does  reveal  that  the  stack  tem¬ 
perature  begins  to  stray  at  the  high  end  of  the  current  draw  profile, 
suggesting  that  the  controller  has  difficulty  in  properly  adjusting 
the  cooling  water  flow  rate.  Considering  the  controller  model  and 
the  correlation  between  the  current  draw  and  necessary  cooling 
water  flow  rate,  it  would  be  expected  that  the  controller  perfor¬ 
mance  would  degrade  when  the  second  and  third  terms  were  too 
small.  If  the  adjustments  made  to  the  flow  rate  were  too  small,  then 
the  stack  temperature  would  be  expected  to  rise,  which  is  exactly 
the  result  seen.  This  breakdown  begins  to  occur  at  a  current  draw  of 
about  280  A,  which  is  beyond  the  applicability  of  the  polarization 
curve  correlation. 


Fig.  3.  Stack  dynamics  for  linear  current  draw. 


Fig.  4.  Cooling  water  dynamics  for  linear  current  draw. 
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Fig.  5.  Stack  dynamics  for  pulse  function  current  draw. 
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Fig.  6.  Cooling  water  dynamics  for  pulse  function  current  draw. 


Fig.  7.  Stack  dynamics  for  high  current  draw. 
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Fig.  8.  Cooling  water  dynamics  for  high  current  draw. 


Fig.  9.  Stack  dynamics  for  low  current  draw. 


The  next  simple  current  profile  considered  is  that  of  a  pulse 
function.  In  this  case,  it  is  possible  to  examine  the  ability  of  the 
controller  to  react  to  large,  sudden  changes  in  the  current  draw. 
Sudden  changes  such  as  this  are  fairly  common  in  many  real  world 
applications,  so  it  is  vital  that  the  controller  be  able  to  quickly  and 
effectively  respond.  The  results  of  this  model  are  given  in  Figs.  5  and 
6,  and  demonstrate  that  the  designed  controller  is  very  effective  at 


handling  isolated  step  changes.  The  stack  temperature  shows  a 
visible  controller  response  within  about  50  s  of  each  disturbance, 
allowing  the  stack  temperature  to  change  by  only  about  2.5  I< 
before  returning  to  the  target  temperature  in  each  case.  That  cor¬ 
responds  to  a  97%  reduction  in  the  change  in  stack  temperature 
from  the  reference  case  of  constant  cooling  water  flow.  Fig.  6  also 
provides  the  best  example  of  the  PI  controller  at  work.  At  t  =  4000  s, 
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Fig.  10.  Cooling  water  dynamics  for  low  current  draw. 
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there  is  a  step  change  in  the  current  draw.  This  step  change  causes 
the  controller  to  suddenly  change  the  cooling  water  flow  rate  for 
two  reasons;  the  first  reason  is  simply  the  difference  in  the  steady- 
state  cooling  water  flow,  which  is  the  first  term  in  Eq.  (26).  The 
controller  also  accounts  for  the  deviation  in  the  stack  temperature 
from  the  set  point,  which  explains  the  small  changes  in  the  cooling 
water  flow  rate  for  about  500  s  following  the  step  change  in  current 
draw. 

The  model  was  tested  using  two  current  profiles  obtained  from 
experimental  data  from  a  3  bedroom  house  [4].  The  experimental 
data  was  broken  into  two  current  profiles,  each  with  a  different 
trend.  The  first  profile  considers  the  current  draw  on  a  home  during 
peak  hours,  so  the  profile  has  multiple  drastic  changes  in  the  power 
demand.  The  second  profile  considers  the  power  demand  in  non¬ 
peak  hours,  featuring  sections  of  time  with  little  change  in  the 
current  draw.  The  results  for  the  first  current  profile  based  on 
experimental  data  are  given  in  Figs.  7  and  8. 

The  results  show  that  the  controller  greatly  reduces  the  resulting 
changes  in  stack  temperature.  The  greatest  deviation  in  stack  tem¬ 
perature  is  only  about  1.5  K;  by  comparison,  the  greatest  deviation  in 
stack  temperature  from  the  average  temperature  in  the  model 
without  the  controller  is  approximately  41  K.  That  means  if  the 
model  without  the  controller  was  designed  to  operate  at  is  average 
stack  temperature  (368  K),  then  the  maximum  deviation  from  the 
ideal  temperature  would  be  about  41  K.  Comparing  the  controller 
deviation  in  stack  temperature  to  this  best  case  scenario  model 
without  a  controller  gives  a  96.5%  reduction  in  the  variation  of  stack 
temperature.  The  average  outlet  temperature  of  the  cooling  water  is 
around  320  K,  with  a  maximum  outlet  temperature  of  333  K.  This 
maximum  change  in  the  cooling  water  temperature  could  prove 
problematic  for  a  system  with  recycled  cooling  water,  but  not  for  a 
system  where  the  water  is  disposed  or  used  for  other  purposes. 

The  second  current  profile  from  experimental  data  is  designed 
using  the  same  procedure  as  the  first.  Figs.  9  and  10  show  the 
current  profile  and  the  results  from  the  calculations.  The  results 
show  that  the  controller  is  notably  less  effective,  mainly  because 
the  average  current  draw  for  the  second  profile  is  much  lower  than 
20  A.  The  maximum  deviation  from  350  I<  is  just  under  5  I<  and  the 
deviations  average  about  2  K.  Using  the  same  comparison  as  in  the 
previous  current  draw  profile,  the  average  operating  temperature 
for  the  reference  model  without  the  controller  is  approximately 
335  K,  and  the  maximum  deviation  from  this  temperature  is  48  K. 
This  means  that  the  controller  still  reduces  the  deviation  in  stack 
temperature  by  as  much  as  96%.  Although  the  controller  is  certainly 
less  effective  at  lower  current  draws,  it  still  achieves  the  desired 
result  within  the  given  restraints.  The  required  cooling  flow  rates 
are  not  excessive  in  any  of  the  simulations  and  it  is  not  necessary  to 
use  chilled  water. 

4.  Conclusions 

A  dynamic  model  of  a  water  cooled  fuel  cell  stack  was  developed 
and  simulated  in  conjunction  with  a  state  space  model  of  a  PI 
controller.  The  stack  temperature  dynamics  and  cooling  water 
temperature  dynamics  models  were  linked  and  solved  simulta¬ 
neously,  while  all  other  dynamics  inside  the  fuel  cell  were  consid¬ 
ered  to  operate  at  pseudo-steady-state.  A  time-varying  PI  controller 
was  designed  for  the  theoretical  model  of  a  PEM  fuel  cell  stack 
jacketed  by  cooling  water.  The  controller  adjusts  the  cooling  water 
flow  rate  through  the  jacket  to  maintain  the  stack  temperature  at  an 
optimal  operating  point.  The  model  was  tested  using  multiple  po¬ 
wer  profiles,  including  two  power  profiles  obtained  experimentally 
from  a  residential  house.  It  was  determined  that  the  controller  is 
capable  of  maintaining  the  stack  temperature  within  2  K  for  most 
current  draws,  and  within  5  I<  for  current  draws  below  20  A. 
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Notation 


English  letters 

Aceii:  active  fuel  cell  area,  m2 

Acs',  reactor  cross  sectional  area,  m3 

Q:  lumped  parameter,  W  K1  A-1 

C2:  lumped  parameter,  W  A-1 

Cp:  specific  heat,  J  mol-1  K1 

CV:  conversion  factor 

Dp:  catalyst  particle  diameter,  m 

fs:  saturation  calculation 

F:  Faraday’s  constant,  C  mol-1 

G:  superficial  mass  velocity,  kg  m-2  s-1 

h:  cooling  water  channel  height,  m 

H:  molar  enthalpy,  J  mol-1 

AH298:  standard  heat  of  reaction,  kj  mol-1 

AHrxn:  molar  heat  of  reaction,  J  mol-1 

i:  current  density,  A  m-2 

I:  current,  A 

kc:  controller  proportional  constant,  kg  s-1  K-1 
I:  fuel  cell  stack  length,  m 
m:  mass,  kg 

m:  mass  flow  rate,  kg  s-1 

nceu:  number  of  cells  in  stack 

ne  :  number  of  electrons 

N:  molar  flow  rate,  fuel  cell  stack,  mol  s-1 

p:  partial  pressure,  kPa 

P:  pressure,  kPa 

q:  molar  flow  rate,  reactor,  mol  s-1 

Q;  heat  flow  into  reactor  per  unit  volume,  W  m-3 

Q.:  heat  flow,  W 

r:  reaction  rate,  mol  s-1  m-3 

R:  universal  gas  constant,  J  mol-1  K-1 

T:  temperature,  K 

U:  internal  energy,  J 

UA:  heat  transfer  parameter,  W  K-1 

V:  reactor  volume,  m3 

Vcdi:  cell  voltage,  V 

Vn:  volume  per  cooling  water  zone,  m3 
W:  fuel  cell  stack  width,  m 
W:  work  flow,  W 

X:  fractional  conversion  of  primary  reactant 
Xs:  fractional  excess  oxygen  fed  to  fuel  cell  stack 
y:  gaseous  mole  fraction 

Greek  letters 


e:  efficiency  factor 
p:  viscosity,  kg  m-1  s-1 
p;  stoichiometric  ratio 
extent  of  reaction,  mol  s-1 
p:  density,  kg  m-3 


478 


D.  O’Keefe  et  al.  /  Journal  of  Power  Sources  256  (2014)  470-478 


ij:  transfer  function  dimensionless  parameter 
t 2-  transfer  function  dimensionless  parameter 
tj:  controller  integral  constant,  s 
</>:  void  fraction  of  catalyst  in  reformer 

Subscripts 

c:  catalyst 

cw:  cooling  water 


CO2 :  carbon  dioxide 
C2H5OH:  ethanol 
H2:  hydrogen 
H20:  water 
in:  inlet  condition 
out:  outlet  condition 
sp:  set  point 
st:  stack 

T:  total  sum  of  chemical  components  in  reactor 


